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We present calculations of the energy, pair correlation function (PCF), static structure factor 
(SSF), and momentum density (MD) for the one-dimensional electron gas using the quantum Monte 
Carlo method. We are able to resolve peaks in the SSF at even-integer-multiples of the Fermi wave 
vector, which grow as the coupling is increased. Our MD results show an increase in the effective 
Fermi wave vector as the interaction strength is raised in the paramagnetic harmonic wire; this 
appears to be a result of the vanishing difference between the wave functions of the paramagnetic 
and ferromagnetic systems. We have extracted the Luttinger liquid exponent from our MDs by fitting 
to data around &f , finding good agreement between the exponent of the ferromagnetic infinitely-thin 
wire and the ferromagnetic harmonic wire. 

PACS numbers: 71.10.Hf, 71.10.Pm, 73.63. Nm 
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I. INTRODUCTION 

Landau's theory of Fermi liquids has proven tremen- 
dously successful at describing a wide range of systems 
of interacting fcrmions. In particular, the theory le- 
gitimizes the free electron model by casting fermionic 
systems in terms of weakly-interacting quasiparticles. 
Systems of electrons in ID provide an intriguing ex- 
ample of departure from the Landau Fermi liquid 
paradigm, exhibiting non-Fermi-liquid behavior for any 
finite strength of the electron-electron interaction^ Per- 
haps the simplest model of electrons in one dimen- 
sion is the one-dimensional (ID) homogeneous electron 
gas (HEG), which comprises electrons on a uniform 
positively-charged background^ 

The strong correlation occurring in ID ensures that the 
excitations are not electron-like quasiparticles, but are in- 
stead collective in nature. An appropriate description of 
the low-energy spectrum of the ID HEG comes from the 
theory of Tomonaga and Luttinger j2t— There are several 
experimental signatures of the Tomonaga-Luttingcr (TL) 
liquid which distinguish it from the normal Fermi liquid; 
these are largely accessible to transport and tunneling ex- 
periments. For example, the conductivity of a ID chan- 
nel as a function of the temperature is expected to vary 
logarithmically in the presence of weak disorder for the 
Fermi liquid, and as a power law for the TL liquid^^ 
Analogous relations hold for the differential conductivity 
and the optical conductivity. Also associated with the 
lack of quasiparticles in the TL liquid is spin-charge sep- 
aration, whereby spin and charge excitations propagate 
at different characteristic velocities'^ 

One-dimensional models are easy to envisage, but 
experimental observation of ID behavior is potentially 
problematic. Low-dimensional systems are never en- 
tirely independent of their 3D environment, leading to 
effects which have the potential to obscure the ID be- 
havior. Furthermore, the presence of impurities has been 



shown to alter drastically the behavior of a TL liquid*^ 
However, even in manifestly 3D systems, behavior un- 
ambiguously characteristic of electrons in ID arises sur- 
prisingly frequently. Features associated with the Lut- 
tinger model have been observed in organic conductors 
(e.g., tetrathiafulvalene-tetracyanoquinodimcthane and 
the Bechgaard salts) transition metal oxideS ) 16 i 17 
carbon nanotubes^^ - — edge states in quantum Hall 
liquids j22r— semiconductor heterostructures^ - — con- 
fined atomic gases ^2r— and atomic nanowiresi^ 3 . The- 
oretical work on electrons in ID thus has a large region 
of potential applicability. 

The exactly-solvable Luttinger model describes elec- 
trons moving in one dimension with short-range interac- 
tions and linear dispersion. Studies with long-range in- 
teractions have found that the exponents and excitation 
velocities are nontrivially altered^ One thus expects to 
be able to describe the ID HEG within the Luttinger 
model framework, but the exact behavior of the param- 
eters of the model is largely unclear. The interactions 
that we study here are long-ranged, possessing a 
Coulombic tail. This is most applicable to systems where 
screening is a small effect, such as isolated metallic car- 
bon nanotubes and semiconductor structures where there 
is negligible coupling to the substrate. 

The ID HEG has been studied with a variety of theo- 
retical and computational approaches. The principal dis- 
tinction between various studies is the choice of electron- 
electron interaction. The bare Coulomb interaction, 
1/|jc|, which describes an infinitely-thin wire, is perhaps 
conceptually the simplest choice, although it is largely 
avoided in the literature^ in its original form due to the 
divergence at x = 0. Instead, many previous authors 
have removed the singularity while retaining the long- 
range behavior by investigating interaction potentials of 
the form V(x) oc (x 2 + g? 2 ) -1 / 2 , where d is a parameter 
related to the width of the wire. This interaction has 
been studied analyticall y 34 i 36 and numerically^! 
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Otherwise, one can derive an effective ID interaction 
by factorizing the wave function into longitudinal and 
transverse parts and assuming that the transverse com- 
ponent is the (2D) single-particle ground state of the 
confining potential. The ID interaction is then the ma- 
trix element of the 3D Coulomb interaction with respect 
to the transverse eigcnfunctions i 38 i 39 An example of this 
is the harmonic wire, in which the transverse confine- 
ment is provided by a parabolic potential, leading to a 
Gaussian density profile in the transverse plane. The 
harmonic wire has been studied with quantum Monte 
Carlo (QMC)^ii variants of the Singwi-Tosi-Land- 
Sjolander approach^2r— and the Fermi hypernetted- 
chain approximation.— 

We have studied both the infinitely-thin wire and the 
harmonic wire using QMC. In this article we report 
QMC calculations of the momentum density (MD), en- 
ergy, pair-correlation function (PCF), and static struc- 
ture factor (SSF) of the infinitely-thin wire at a variety 
of densities and system sizes. The MD results in partic- 
ular show the non-Fermi-liquid character of the system 
and allow us to recover one of the TL parameters. The 
total energy data that we provide are exact and may 
be regarded as a benchmark for future work. We also 
present calculations of the MD for the harmonic wire, 
again extracting one of the TL parameters. 

The rest of this paper is structured as follows: the mod- 
els for which we perform our calculations are described 
in Secini In Sec. IIIII we outline QMC methods and pro- 
vide the details of our approach. We report the ground 
state energies of both models in Sec. IIV Al and describe 
the PCFs in Sec. HVBl In Sec. ITVCl we give the SSFs 
that we find for the infinitely-thin wire and in Sec. IIVDI 
we give the MDs for both models. We describe the pro- 
cedure for estimating a parameter of the TL model in 
Sec. IIV El Finally, we draw our conclusions in Sec. |V] 
We use Hartree atomic units (h = \e\ = m e = 4-7T£o = 1) 
throughout this article. 



B. Infinitely-thin wire 

The Ewald interaction for the infinitely-thin wire may 
be written 



n— — OQ 



1 



\Xij + nh\ 



L/2 



dy 



L J-l/2 +nL-y\ 



(2) 



which is calculated in practice using an accurate approx- 
imation based on the Euler-Maclaurin summation for- 
mula; see Eq. (4.8) of Ref. E3 for details. 

The interaction of Eq. ([2]) diverges as when 
Xij — > 0. In higher dimensions, the divergence in the 
interaction energy is canceled by an equal and opposite 
divergence in the kinetic energy, so that nodes do not 
necessarily occur where two antiparallel spins occupy the 
same position^ In the infinitely-thin ID system, the cur- 
vature of the wave function is unable to compensate for 
the divergence in the interaction potential, so the trial 
wave function has nodes at all of the coalescence points 
for both parallel and antiparallel spin pairs. The result is 
that the ground state energy is independent of the spin- 
polarization and depends only on the density. In other 
words, the Licb-Mattis theorem^ 9 - does not apply and the 
paramagnetic and ferromagnetic states are degenerate for 
the interaction of Eq. @. We have examined only the 
fully spin-polarized case for the infinitely-thin wire. 



C. Harmonic wire 

The second model we have studied describes electrons 
in a 2D confinement potential given by 



V±(r x ) =ri/86 4 



(3) 



II. MODELS 



A. Hamiltonian 



The Hamiltonians for both of the models we have stud- 
ied may be written as 



H 



1 N rl 2 



i=i 



N 

) + TT^Mad , (1) 



where V^iad is the Madclung energy (the interaction of a 
particle with its own background and periodic images), 
X -ij — | X £ Xj | IS the distance between electron i and 
electron j, and V(xij) is the Ewald interaction; this is 
the interaction of an electron at Xi with another electron 
at Xj, all of electron j's periodic images, and 1/iV-th of 
the uniform positive background. The two models that 
we have studied differ in the V{xij) and VmoA terms. 



where b is the width parameter and r is the magnitude 
of the projection of the electron position onto the plane 
perpendicular to the axis of the wire. The Ewald-likc 
interaction for this model may be written a o 39 ' 50 
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\xij - mL\ 



erf , \ x ij - mL \ 



2b 



Ei [{bGnf] cos(G 



(4) 



71=1 



where G = 2-k/L. Equation (|4| possesses a long-range 
Coulomb tail and is finite at x^ = 0. A derivation of Eq. 
(|4"| is given in Appendix [XJ For the harmonic wire we 
have probed different polarizations, £ = — N±\/N. 
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III. DETAILS OF CALCULATIONS 

We computed expectation values using the varia- 
tional and diffusion Monte Carlo (VMC and DMC, 
respectively) methods as implemented in the CASINO 
program^ For the infinitely-thin wire, we combined 
VMC and DMC results to form extrapolated estimates 5 ^ 
where applicable, whereas for the harmonic wire we used 
VMC alone. 

In the VMC method the expectation value of the 
Hamiltonian with respect to a trial wave function is cal- 
culated using a stochastic integration technique.— Trial 
wave functions usually contain a number of free param- 
eters; we optimized the free parameters in our wave 
function by unreweighted variance minimization 5 ^— and 
linear-least-squares energy minimization^ DMC is a 
stochastic projector technique for solving the many-body 
Schrodinger equation and generates configurations dis- 
tributed according to the product of the trial wave func- 
tion and its ground state component . 52 ' 57 DMC calcu- 
lations of expectation values of operators that commute 
with the Hamiltonian are in principle exact for systems 
in which the wave function nodes are known; this is the 
case for both the infinitely-thin and harmonic wires. 

We used a Slater-Jastrow wave function for both sys- 
tems, where the Jastrow factor comprised two-body 
terms consisting of smoothly truncated polynomials and 
a sum of cosines with periodicity commensurate with that 
of the simulation cell<££ The orbitals in the Slater deter- 
minant were plane waves with wave vectors up to fcp = 
7r/(4r s ) for the paramagnetic systems and kp = ■k/{2t s ) 
for the ferromagnetic systems. The orbitals were evalu- 
ated at quasiparticle coordinates related to the actual co- 
ordinates by a backflow transformation.— Backflow pro- 
vides an efficient way of describing three-body correla- 
tions in the ID HEG, but leaves the exact nodal surface 
unchanged. 

One method for assessing the wave function quality 
is to examine the fraction of the correlation energy re- 
trieved, (£hf - £vmc)/(£hf - Sdmc), where £hf is 
the Hartree-Fock energy, and E'dmc and -Evmc are the 
DMC and VMC energies, respectively. We tested several 
types of wave function for the infinitely-thin wire with 
r s = 15 a.u., N = 15, and £ = 1; our VMC calculations 
retrieved 99.9989(9)% of the correlation energy when we 
used a two-body Jastrow factor and backflow transforma- 
tions [the error bars were O(10~ 8 ) a.u.], which is the type 
of wave function we use throughout this paper. While it 
is indeed the case that DMC is formally exact for the ID 
HEG, the quality of the trial wave function is important 
for the statistical efficiency of the DMC method and the 
accuracy with which expectation values of operators that 
do not commute with the Hamiltonian may be computed. 



The DMC energy did not change beyond statistical error 
upon varying the number of walkers between 640 and 
2000, so we used ~ 1000 walkers in our calculations and 
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FIG. 1: DMC energy of the infinitely-thin wire for several 
different timesteps and a linear fit to the data. The plot is for 
r 3 = 1 a.u. and N = 37 with 1000 configurations. 



assumed that population control bias is negligible. The 
dependence of the energy upon the DMC timestep r was 
also investigated; Fig.Q]sIiows that for small r the energy 
is constant. We performed our calculations at a single 
timestep given by r = 0.008 r 2 . This fairly conservative 
choice was made to ensure that time step bias is entirely 
negligible. The RMS distance diffused by each electron 
in a single step was thus slightly less than r s /10. 

For the infinitely-thin wire, we used simulation cells 
containing 37, 55, 73, and 99 particles subject to periodic 
boundary conditions for our calculations of the energy, 
PCF, and SSF. Our MD calculations for the infinitely- 
thin wire also used a much larger cell with N = 255. 
For the harmonic wire, we used cells with ./V = 123, 155, 
and occasionally 255 for the £ — 1 systems and cells with 
N = 22 and 102 for the C = systems. 

Previous work encountered difficulties in sampling dif- 
ferent spin configurations of the harmonic wire for ( =/= 1 
due to the presence of "pseudonodes" at the antiparal- 
lel coalescence points^ although these problems were 
largely overcome by the use of lattice-regularized diffu- 
sion Monte Carlo (LRDMC) in Ref. |H The problem 
occurred because for strong, repulsive interactions the 
wave function can become small when two antiparallcl 
spins approach one another. Combined with a small 
time step this can lead to simulations where opposite 
spins exchange positions infrequently and the space of 
spin configurations is explored very inefficiently. Use of 
a small time step is a necessary part of the algorithm of 
projector methods like DMC. We have avoided crgodic- 
ity problems by using VMC to study the harmonic wire; 
in the VMC method there is no restriction other than 
ergodicity on the transition probability density and one 
may propose moves however one wishes provided that the 
acceptance probability is modified accordingly. We use 
electron-by-electron sampling with the transition proba- 
bility density given by a Gaussian centered on the initial 
electron position. The VMC "time step" in fact bears 
no relation to real time and is simply the variance of the 
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r s (a.u.) 


b (a.u.) 


^exch 


1 


0.1 


0.051(1) 


1 


4 


0.160(2) 


15 


0.1 


0.0016(2) 


15 


4 


0.0020(3) 



r s (a.u.) N £dmc (a.u. / elec. 



TABLE I: Frequency with which electrons' paths cross in our 
VMC simulations of the harmonic wire. The quantity s cxc h is 
the proportion of proposed single-electron moves that result 
in a change in the ordering of the particles. A typical calcula- 
tion comprises between 10 7 and 10 s proposed single-electron 
moves. The data shown are for N = 22. 
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TABLE II: The DMC energies for the infinitely-thin wire. 



FIG. 2: (Color online) Plot of the DMC energy against the 
reciprocal of the square of the system size for the infinitely- 
thin wire. The energy has been offset by the extrapolate E x 
that one obtains using the form E(N) = Eac + BN~ 2 . 

transition probability density. In practice, the unmodi- 
fied time steps (chosen to achieve a 50% acceptance ratio) 
used in VMC are usually large enough to eliminate er- 
godicity problems, although wc found some cases where 
it was necessary to enforce a lower limit on the width 
of the transition probability density. Table Q] shows the 
frequency with which electrons changed positions in our 
simulations for both high and low density systems with 
strong and weak confinement. 

IV. RESULTS 
A. Energies 



For the infinitely-thin wire, we used DMC to calculate 
the exact ground state energy since there is no ergodicity 
problem. Tabic [XT] shows the DMC energies that we ob- 
tained for r s = 1, 2, 5, 10, 15, and 20 a.u., with N = 37, 



r s (a.u.) Eoo (a.u. / elec.) 

1 0.1541886(2) 

2 -0.20620084(7) 
5 -0.20393235(2) 
10 -0.142869097(9) 
15 -0.110466761(4) 
20 -0.090777768(2) 



TABLE III: The DMC energies for the infinitely-thin wire 
extrapolated to the thermodynamic limit using the form 
E(N) = Eoo + BN~ 2 . 



55, 73, and 99 particles. We find that the Fourier trans- 
form of the two-body Jastrow factor u(k) takes the form 
u(k) oc 1/k as k — > 0, allowing us to estimate the leading- 
order scaling of the finite-size correction to the kinetic 
cncrgyj££ Furthermore, we observe that the static struc- 
ture factor S{k) goes to zero linearly at k = 0, allowing 
calculation of the corresponding correction for the po- 
tential energy^ Motivated by these results, we use the 
form E(N) = + BN~ 2 to extrapolate to the ther- 
modynamic limit. Figure [5] demonstrates that this form 
fits the data well, and Table IIIII shows the extrapolated 
energies E^. 
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0.0901489(7) 


0.09014(1) 
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-0.1631207(8) 


-0.16311(2) 


10 


-0.1231560(3) 


-0.123157(3) 


15 


-0.0971194(1) 


-0.097120(2) 


20 


-0.0807160(2) 


-0.080717(1) 



TABLE IV: Comparison of our VMC energies for the har- 
monic wire (6=1 a.u., f = 1) with those of Ref. |39| . calcu- 
lated using the LRDMC method. For both sets of results the 
energies were extrapolated to the thermodynamic limit using 
the functional form E(N) = E^ + BN' 1 + CN~ 2 , where E x , 
B, and C are fitting parameters. 



The many-body Bloch theorem states that the wave 
function ipx satisfies^ 



^t(xi,. . .,Xj+L, . 



,xn) , (5) 



where L is the length of the simulation cell and k a is the 
simulation cell Bloch wave number. Averaging over k a 
in the irreducible Brillouin zone ("twist averaging") has 
been shown to reduce greatly single-particle finite-size ef- 
fects in two and three dimensions i£2r— In ID, however, 
use of a nonzero k a does not result in reoccupation of the 
orbitals, merely adding k a /2 to the energy per particle 
and leaving unchanged (or trivially altering) other expec- 
tation values. It is easy to show that the average of k a /2 
over the Brillouin zone falls off as 0(N~ 2 ) in ID; hence 
single-particle finite-size effects simply lead to an addi- 
tional 0(N~ 2 ) error in the energy per particle, which is 
removed when we extrapolate to infinite system size. 

The trial wave functions in our calculations for the har- 
monic wire are of sufficient quality that the variational 
energies we obtain are in statistical agreement with exact 
results in the literature^ Table IIVI shows the compari- 
son. 
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FIG. 3: (Color online) PCF of the infinitely-thin wire at sev- 
eral densities. The data shown are for N = 99 and are ex- 
trapolated estimates \2,guyic{x) — gvMc(x)]. 
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FIG. 4: (Color online) PCF of the harmonic wire with r s — 1 
a.u., N — 39, b = 1 a.u., and £ = 1. The solid line shows 
our VMC results and the symbols show the LRDMC results 
of Ref. [3!| The inset shows the same data at the origin. 



configurations generated by the QMC algorithms. The 
antiparallel-spin PCF may be written as 



The PCF is accumulated in QMC simply by binning the 
interparticle distances throughout the simulation. The 
parallel-spin PCF is 

1 / N ' \ 
SttO) = — \ ^2 s (\ x i,t - x iM ~ x )) . ( 6 ) 
\i>j I 

where p a is the average density of electrons with spin 
tr, Xi a is the position of the ith electron with spin a 
and the angular brackets denote an average over the 



9n( x ) = (y2y2 S (\ x i,t- x u\- x )) ■ ( 7 ) 

The PCF for the harmonic wire was calculated for dif- 
ferent confinements and system sizes by Casula et al. us- 
ing the lattice-regularized DMC method.— Figures|4]and 
E] show the agreement of the LRDMC results with the 
present work. Figure [3] shows the PCF for the infinitely- 
thin wire at several values of r s . 
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FIG. 5: (Color online) PCF of the harmonic wire with r s = 1 
a.u., N = 42, b — 1 a.u., and £ = 0. The solid line shows 
our VMC results and the symbols show the LRDMC results 
of Ref . [39l. The function plotted is [g^(x) + g^i(x)]/2. 
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FIG. 6: (Color online) SSF of the infinitely-thin wire at sev- 
eral system sizes. The data shown are extrapolated estimates 
[2SbMc(fc) — SvMc(k)] for r s — 2 a.u.. The main plot shows 
the behavior at the peak and the inset shows a zoomed-out 
view. The PCF was not extended beyond L/2. 



FIG. 7: (Color online) Effect of extending the PCF before 
performing the transformation of Eq. JjjJ. The square symbols 
(labeled "WS radius") show the SSF obtained from the finite- 
cell PCF. The solid, dash-dot, and dashed lines (which all lie 
on top of one another) are from the N = 99, 73, and 55 PCFs, 
respectively, where in each case the PCF has been extended 
out to many simulation cell lengths using the fitting form of 
Eq. ((9j). The data shown are for r 3 = 20 a.u. 
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FIG. 8: (Color online) SSFs of infinitely thin wires obtained 
from Eq. JHJ and the extended N = 99 PCFs. 



Static structure factor 



The SSF of the ID HEG is defined asi 

S(k) = l + ^[ \g{x) - l]e- lk * dx , (8) 



and SSFs that we present here are for the ferromagnetic 
infinitely-thin wire. As explained in the introduction, the 
antiferromagnctic and ferromagnetic phases are degener- 
ate for the infinitely-thin wire, so we do not violate the 
Lieb-Mattis theorem with our choice of system. 

The PCF can only be directly measured in QMC for 
x < L/2 due to the finite extent of the simulation cell. 
This manifests itself as a scaling of the SSF peak at 2fcp 
with system size, as demonstrated by the data in Fig. [6] 
The height of the 2fcp peak in the finite-cell SSFs does not 
scale as N (and so L) to any single power but appears to 
be sub-linear, consistent with the presence of quasi long- 
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FIG. 9: (Color online) SSF of the infinitely-thin wire at very 
low densities. The results shown are for N = 99. The finite- 
cell VMC PCFs were used to generate the SSFs in the plot. 



range order. At k away from the peak the SSFs appear 
to agree very well for different cell sizes. 

We further investigated finite-size effects by perform- 
ing a fit to the oscillatory tails of the PCF and using the 
fitted function to extend the PCF far beyond L/2 be- 
fore using Eq. (j8]) to calculate the SSF. After testing a 
number of functional forms, we found that a good-quality 
and simple fit to the oscillatory tails of the PCF takes the 
for m 34 i 39 



g(x) — 1 = Acos(2fcpa;) exp(- 



(9) 



where A and B are treated as fitting parameters. The 
choice of Eq. © is motivated by the charge-charge cor- 
relation function of Ref . [34l The parameters we obtained 
when fitting Eq. (J9)) to our results are given in Table [V] 
in Appendix iBl 

We fitted Eq. © to the PCF data Idi- 
ot's < x < L/2 — 6r s , although wc found that the 
results were not very sensitive to the region of data 
included in the fit. The data close to the origin were 
not included in the fit since Eq. © is only a good fit 
for long-range correlations. The data at the edge of the 
cell were excluded because that is the region midway 
between the electron at the origin and its next periodic 
image, and might be expected to be a region where the 
PCF suffers particularly badly from finite-size effects. 

We then formed the extended PCF by reinstating all 
of the original PCF data up to L/2 — 6r s and append- 
ing a tail for x > L/2 — 6r s using Eq. ^ and the fitted 
parameters. Performing the Fourier transform of Eq. J8J 
numerically on the extended PCF results in a SSF (fai- 
rs < 20 a.u.) with a greatly-enhanced peak at 2/cp, but 
that agrees very well with the finite-cell SSF everywhere 
else. Figure [7] shows the difference between the SSFs ob- 
tained from the finite-cell and the extended PCFs. Under 
the extension scheme, the peak at 2fcp appears to be sus- 
ceptible to noise; in particular, the density of fc-points 



FIG. 10: (Color online) MD of the infinitely-thin wire at sev- 
eral densities. The data shown are for N = 99 and are extrap- 
olated estimates [2nuMc(k) — nvMc(&)]- The statistical error 
bars are much smaller than the symbols and some symbols 
have been omitted for clarity. 



at which the SSF is calculated heavily affects the appar- 
ent height. The fitting function of Eq. © possesses a 
peak at 2kp in Fourier space and smoothly decays away 
elsewhere, and is expected to be closer to the N = oo 
limit. 

The asymptotically-correct charge-charge correlation 
function of Schula 3 -! includes higher order terms contain- 
ing oscillations at wave numbers given by even multiples 
of fcp- For r s < 15 a.u. we find no discernable features 
at larger k. However, a small feature at 4fcp starts to 
develop at r s sa 15 a.u., and for r s = 20 a.u. we ob- 
serve a clear peak, visible in Figs. [JJ [51 and [HI We per- 
formed short VMC calculations at extremely low den- 
sities, r s = 50 a.u. and 100 a.u., where the electron- 
electron coupling is very large, to search for more notice- 
able features at k > 2kp. We find that peaks in the SSF 
do indeed appear at even multiples of kp for these sys- 
tems, as evidenced in Fig. The SSF of the r s = 100 
a.u. system has clear peaks at k = 2fcp, 4/cp, and 6fcp. 
This suggests that one could add higher order terms to 
the fit of Eq. © for the low density systems and perform 
the extension scheme again, although this seems unlikely 
to produce any new behavior. 



D. Momentum density 



The MD is accumulated as 



n(k) = (- 



* T (r) 
*t(xi) 



exp[z/c(xi — r)] dr 



(10) 



8 




FIG. 11: (Color online) VMC MD of the harmonic wire with 
b = 0.1 a.u. and £ = at several densities. The data shown for 
each density are for N = 22 and 102 (joined to form one data 
set). The statistical error bars are smaller than the symbols. 




0.2 
A: (a.u.) 



will be singular^ 

For the systems with £ = 1, we have used kp = 
7r/(2r s ), whereas for the systems with £ = 0, we 
have used kp = 7r/(4r s ). Figure [TU] shows the MDs 
that we obtain by evaluating the extrapolated estimator 
2^DMc(fc) — ^vmc(^) f° r t ne infinitely-thin wire. The 
VMC and DMC results differed by no more than ~ 2 
error bars, so that evaluating the extrapolated estimator 
changed the results very little. Figure [Til shows the MD 
for the harmonic wire with b = 0.1 a.u. and ( = 0. 

A particularly interesting feature of the paramagnetic 
harmonic wire MD is that as r s is increased and b is de- 
creased the function shifts much of its weight to larger 
k, and n(0) reduces to values around 0.5. This is a di- 
rect manifestation of the harmonic wire becoming more 
like the ferromagnetic infinitely-thin system. One can in 
some cases see a feature resembling the gradient discon- 
tinuity appearing at 7r/(2r s ), i.e., at twice the paramag- 
netic Fermi wave vector. In particular, for r s = 10 a.u. 
and b = 0.1 a.u. the MD possesses a feature at ?r/(2r s ). 
Upon closer inspection we find that the MD for the un- 
polarizcd system with b = 0.1 a.u. agrees very well with 
that of the infinitely-thin wire (6 = and £ = 1). Figure 
1121 illustrates this comparison. It thus appears possible 
to in some sense tune the effective Fermi wave vector 
by adjusting the strength of the confinement (and the 
density). A dense paramagnetic system with very weak 
confinement shows significant occupation of momentum 
states up to approximately 7r/(4r s ). Increasing the ef- 
fective interaction strength increases this value of k until 
it eventually saturates at the ferromagnetic kp- This 
reflects the fact that in the limit r s — > oo the pseudon- 
odes at antiparallel-spin coalescence points become true 
nodes. 



E. Tomonaga-Luttinger liquid parameters 



FIG. 12: (Color online) Comparison of the r s — 10 a.u. MD 
for the infinitely-thin wire (6 = 0, £ = 1) with that obtained 
for the harmonic wire with b = 0.1 a.u. and £ = 0. The 
statistical error bars are similar in size to the symbols. The 
dashed lines show the values of k at which Eq. pip was fitted 
to the data for the calculations of the exponent a. 



where 'FtM is the trial wave function evaluated at 
(r, X2, ■ ■ ■ , x n ) and angular brackets denote an average 
over configurations. The MD is the integral of the spec- 
tral function from minus infinity up to the chemical 
potential^ The MD exhibits a drop at k-p because the 
peak in the spectral function reaches and passes through 
the chemical potential. If the peak in the spectral func- 
tion is a (5-function at kp {i.e., the spectral function pos- 
sesses a quasiparticle peak) then the MD is discontinuous 
at the Fermi edge. However, in ID we expect the excita- 
tions to be collective rather than single-particlc-like. The 
ID systems should thus have MDs that are continuous at 
kp, although TL liquid theory predicts that the gradient 



Close to the Fermi wave vector, TL liquid theory suggests 
that the MD takes the forrn^ 

n{k) = n(kp) + A[sign(fc - k F )]\k - kp\ a , (11) 

which we have fitted to our results treating n(kp), A, 
and a as fitting parameters. Note that within TL liq- 
uid theory the exponent a is related to the TL liquid 
parameter— K p by 
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FIG. 13: (Color online) Exponent a [found from fitting Eq. 
(|lf \ to the MD] against the range of data included in the 
fit. The range of data is described by \k — kp\ < £&f. The 
symbols are the fitted exponent values and the solid lines are 
linear fits to the exponents in the region e > 0.05. The data 
shown are for the infinitely-thin wire. 
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FIG. 14: (Color online) Exponent a found from fitting Eq. 
(|If \ to the MD around fc = n/(2r 3 ) for the £ = 1 systems 
and k = 7r/(4r s ) for the £ = systems. 
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FIG. 15: (Color online) Exponent a found from fitting Eq. 
(JTTJ to the MDs of the ferromagnetic systems. 




FIG. 16: (Color online) Exponent a found from fitting Eq. 
(|f f P to the MD of systems with r 3 — 1 a.u. 



If the range of data included in the fit is described by 
| k — fcp | < e&F , the choice of e can present some difficul- 
ties. Ideally, one would choose e —> since Eq. (fTT|) is 
potentially valid for k — > fcp, and indeed using the entire 
range of MD results yields rather poor fits. However, the 
estimate of a becomes noisy when e is small, and at the 
extreme where just two data points are included, one can 
of course obtain any value for a. This leads us to include 
fits constructed using a larger range of k values. In prac- 
tice, we chose to perform a linear extrapolation to e = 
excluding fits where e < 0.05 for the £ = 1 systems, as 
shown in Fig. [T31 For the £ = systems, the extrapola- 
tion used a from fits for which e > 0.25. The trend that 
we observe in the exponent with respect to e is similar 
to that found in Ref. WX 

Figure Q3] shows the exponents a that we obtain for 
several densities, polarizations, and confinements. All 
of the systems show the same general trend; a tends to 
in the high-density limit and to 1 in the low-density 
limit. As mentioned earlier, it is important to note that 
for the C = 1 systems we fitted Eq. (fTTj) to the MD at 
7r/(2r s ), whereas for £ = we used 7r/(4r s ). The change 
in shape of the MD upon varying the interaction strength 
that we noted in Sec. lIVDl and the apparent shift in /cf, 
suggests that one could also extract a relevant exponent 
from fitting to other values of k. For example, we showed 
in Fig.fT^the similarity between the r s = 10 a.u., b = 0.1 
a.u., £ = MD and that of the r s = 10 a.u., 6 = 0, 
C = 1 system. Despite the similarity of the MDs for 
the two systems, the fits used to extract a from each 
were performed at different values of k — a factor of 
two apart in fact. The result is that the exponent for 
the paramagnetic wire is larger, since the Fermi edge for 
that system has apparently shifted to k above the fitting- 
region. 

With this in mind, Fig. [15] shows the £ = 1 results 
alone, since for the ferromagnetic systems one can clearly 
and reliably state that &f = 7r/(2r s ) for the whole range 
of densities. The exponent a for the infinitely-thin wire 
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is reasonably well-approximated by the function 
a = tanh(r s /8) , 



(13) 



which gives a maximum deviation of 0.011(3) from the 
b = QMC results, which occurs at r s = 15 a.u. The 
exponents for the harmonic wire with 6=1 a.u. and £ = 
1 show a maximum deviation from Eq. (|13[) of 0.057(6), 
which we find at r s = 5 a.u. 

The exponent a has been reported in previous theo- 
retical and experimental studies. Reference [5(3 gives the 
exponents for b = 0.1, 1, and 4 a.u. (with r s = 1 a.u. 
and £ = 0) from VMC calculations. In Fig. [1^] we have 
shown how the results given there compare with ours. 
It appears that the principal difference between the two 
studies is the procedure for deciding upon a fitting re- 
gion; Ref. [Ho| docs not give details of any extrapolation 
to e = and presumably the whole range of n(k) was 
included in the fit. Figure [TBI also includes the exponent 
we find for the infinitely-thin wire (from VMC and DMC 
estimates of the MD) at b = 0. 

The exponent a has also been reported from ex- 
periments, mostly through measurements accessing the 
single-particle density of states near the Fermi edge. The 
exponent for carbon nanotubes ranges between 0.2 and 
0.4, although it is difficult to map the behavior of elec- 
trons in these systems onto our model since the elec- 
tronic properties depend on the folding geometry^i^— 
For the Bechgaard salts, which have a ID carrier density 
of r s sa 6.9 a.u., exponents between 0.5625 and 0.8 have 
been report ed. 11 ' 68 " — 



V. CONCLUSIONS 

We have presented calculations of the ground state 
energy, PCF, SSF, and MD of the infinitely-thin ID 
HEG model using VMC and DMC. We observe the 
development of peaks at increasingly-large even-integer- 
multiples of kp in the SSF as the density is lowered, con- 
sistent with the predictions of Schulz<2i 

For the harmonic wire model, we have reported 
ground-state MDs and TL parameters for a range of den- 
sities and confinements. We used VMC to produce these 
results; comparison of our PCFs, SSFs, and ground state 
energies with LRDMC results^ where available indicates 
that our results are extremely accurate. 

The MDs of the £ = systems tend towards the MDs 
of the infinitely-thin wire and ferromagnetic harmonic 
wire as b is decreased and as r s is increased, both of 
which have the effect of increasing the electron-electron 
coupling. One interpretation for this is that correlation is 
dominating over kinetic confinement, so that antiparallel 
spin pairs are avoiding one another almost as much as 
parallel spin pairs. 

The TL parameters calculated for the b = 1 a.u., £ = 1 
system show reasonable agreement with the infinitely- 
thin wire results; the maximum deviation that we observe 
between the parameters for the two systems is 0.051(6), 



which occurs at r s = 5 a.u. The exponent a, which 
describes the behavior of the MD at fcp, takes values 
between and 1. The exponent for the £ = systems 
shows the same general trend, although the value of a 
is typically higher than for the £ = 1 systems. This is 
largely a consequence of the shift of the weight in the MD 
(including the singularity in the gradient) to larger k as 
the coupling is increased. 
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Appendix A: Derivation of the quasi-lD interaction 

One may derive Eq. Q from first principles. Suppose 
the energy scales are such that we may write the wave 
function as a product 9(r±)ip(x), where x is the projec- 
tion of the electron position onto the axis of the wire and 
r*j_ is the transverse position. 

If the electrons are sufficiently confined in the trans- 
verse plane, one may obtain the ID interaction v(x) by 
integrating over the transverse part of the wave function, 



v(x) 



l%±)| 2 |0(r'x)| s 



\r± — r 



I 12 



1/2 



dr\dr'\ 



(Al) 



For the harmonic wire, the confining potential is r^/86 4 , 
where b is a parameter. If r s 3> 7r6/4, one may make 
the assumption that the electrons occupy only the lowest 
sub-band, which is given by 



0(r ± ) = 



1 



: exp 



Substituting Eq. (J52J) into Eq. (TAT]) yields^ 2 - 



(A2) 



(A3) 



which is finite at x = but retains a long-range l/\x\ 
tail. The Fourier transform of Eq. (| A3[) is 



v(k) = #i(& 2 fc 2 )exp(& 2 fc 2 ) 



(A4) 



where E\ is the exponential integral function. 

Having found the real and reciprocal space representa- 
tions of the ID interaction in a harmonic wire, one must 
perform an Ewald-likc sum to enable calculations with 
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periodic systems. We follow a route similar to that of 
Ref. M 

The interaction of an electron at the origin with an- 
other at position x, all of that electron's periodic images, 
and its background is given by 



which is a form suitable for numerical evaluation. 

One may take the first step towards simplifying 72 (x) 
by performing a Poisson summation on p{y) , 



m— — c>D 



)= V \v(x-mL)-- dy v(x-mL-y)\ , where G = 2ir/L. Putting Eq. (|A14|) into Eg. (| A10[) gives 



p{y) = j 



l + 2j2^ {bGn)2 cos(Gny) 



(A14) 



(A5) 

where L is the length of the simulation cell. The objec- 
tive is to write Eq. (|A5[) in terms of quickly converging 
discrete sums. The first step is to write Eq. (|A5|) in the 
more useful form 



0) = j (x) - - 



dy v(x - y) , 



where 



To 



(x) = v(x — mL) 



(A6) 



(A7) 



m— — 00 



Equation (|A7|I is already in a form that is quick and easy 
to evaluate, so we turn our attention to reformulating the 
integral in the second term of Eq. (|A6|) . We first perform 
the trick of both adding and subtracting a Gaussian term 
P(y), giving 



dy v(x — y) = 71 (x) + 72 (x) , 



where 



7i (z) = - 



dy v(x - y)p{y) , 



p{y) - -j- 



(A8) 

(A9) 
(A10) 



72 (x) = / dyv(x-y) 
J — 00 

and the term that we have added and subtracted is 
00 1 / 1 \ 

p(»)= E 2bj^ exp {-4tf {y - mL)2 ) ■ (A11) 

rn— — 00 V \ / 

It is clear that <j>(x) may now be written simply as 

<j)(x) = 70 (x) +7i(x) +72(x) . (A12) 

We first inspect 71 (x), finding that it may be integrated 
directly to give 



00 , 1 

m— — 00 v 



|x — mL\ 
2b 



(A13) 



72 



,-(6Gr 



w(x — y) cos(Gny) , 



(A15) 

which may straightforwardly be rewritten in its final 
form, 



72 (x) = 



2V2^ 
L 



v{Gn) cos(Gnx) 



-(bGnr 



(A16) 



n=l 



where we have used the result 



c?j/ v(x — y) cos(Gny) = y 2tt v(Gn) cos(Gnx) . 

(A17) 

Finally, putting the expressions for the 7 functions, 
Eqs. (|A"7|) . (|A13|) . and (|A16|) . into Eq. (|A12|) and remem- 
bering that is given by Eq. (|A4[) . we obtain the more 
computationally convenient form 



X — 77l_L| 

26 



1 



|x — mXI 



erf 



x — mi I 
26 



I 



^ Si [(6Gn) 2 ] cos(Gnx) 



(A18) 



n=l 



It should be noted that in Ref. |39|, Rydberg rather than 
Hartree units were used so that the potentials given there 
differ from Eq. (|A18[) by a factor of 2. 



Appendix B: Pair correlation function fitting 
parameters 

Table [V] shows the parameters that we obtained when 
fitting Eq. © to the extrapolated estimates of the PCF 
for the infinitely-thin wire. We performed the fit for r s = 
1, 2, 5, 10, 15, and 20 a.u. with systems containing = 
37, 55, 73, and 99 particles. The PCF data in the range 
6r s < x < L/2 — 6r s were included in the fit. 
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r s (a.u.) 


N 


A 


B (a.u.) 


16 


1 


37 


1.908 


3.291 




1 


55 


3.090 


3.683 












17 


1 


7:5 


1.940 


3.446 




1 


99 


2.113 


3.671 




2 


37 


4.851 


2.952 


18 


2 


55 


4.573 


2.979 




2 


73 


6.047 


3.069 




2 


99 
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3.273 




5 


37 
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1!) 


5 


55 
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5 


73 


10.061 
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20 


5 


99 
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10 


37 
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10 


55 
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21 


10 


73 
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1.839 


22 


15 


37 
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25 
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37 
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20 


55 


8.361 
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20 


73 
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2(1 


20 


99 


8.895 


1.396 





TABLE V: Table showing the fitting parameters A and B 
from Eq. ([9]) obtained from fitting to the extrapolated esti- 
mates of the PCF for the infinitely-thin wire. The fits were 
to PCF data in the range 6r s < x < L/2 — 6r a . 
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